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LIOUVILLE THEOREMS AND SPECTRAL EDGE 
BEHAVIOR ON ABELIAN COVERINGS OF COMPACT 

MANIFOLDS 

PETER KUCHMENT AND YEHUDA PINCHOVER 


Abstract. The paper describes relations between Liouville type 
theorems for solutions of a periodic elliptic equation (or a system) 
on an abelian cover of a compact Riemannian manifold and the 
structure of the dispersion relation for this equation at the edges 
of the spectrum. Here one says that the Liouville theorem holds 
if the space of solutions of any given polynomial growth is hnite 
dimensional. The necessary and sufficient condition for a Liouville 
type theorem to hold is that the real Fermi surface of the elliptic 
operator consists of finitely many points (modulo the reciprocal 
lattice). Thus, such a theorem generically is expected to hold at 
the edges of the spectrum. The precise description of the spaces 
of polynomially growing solutions depends upon a ‘homogenized’ 
constant coefficient operator determined by the analytic structure 
of the dispersion relation. In most cases, simple explicit formulas 
are found for the dimensions of the spaces of polynomially growing 
solutions in terms of the dispersion curves. The role of the base of 
the covering (in particular its dimension) is rather limited, while 
the deck group is of the most importance. 

The results are also established for overdetermined elliptic sys¬ 
tems, which in particular leads to Liouville theorems for polynomi¬ 
ally growing holomorphic functions on abelian coverings of compact 
analytic manifolds. 

Analogous theorems hold for abelian coverings of compact com¬ 
binatorial or quantum graphs. 


1. Introduction 

The classical Liouville theorem claims that any harmonic function 
(i.e., a solution of the Laplace equation Aw = 0) in M” that has a 
polynomial upper bound is in fact a (harmonic) polynomial. In par¬ 
ticular, the space of all harmonic functions that grow not faster than 
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(7(1 + |x|)'^, is of the finite dimension^ 

no'i h n + N^. /n + iV-2\ 

( 1 - 2 ) hn,N ■ ^ j ^_2 y 

The problem of extending this result to more general elliptic operators 
and to Laplace-Beltrami operators on general Riemannian manifolds 
of nonnegative Ricci curvature^ has gained prominence since the work 
of S. T. Yau |HZj. The questions asked concern the hnite dimensionality 
of the spaces of solutions of a prescribed polynomial growth, estimates 
of (and in rare cases formulas for) their dimensions, and the struc¬ 
ture of these solutions. One can hnd recent advances, reviews, and 
references in 1591 [HU] . In particular, Yau’s conjecture on the valid¬ 
ity of the Liouville theorem for Riemannian manifolds of nonnegative 
Ricci curvature was proven in full generality by T. H. Golding and 
W. P. Minicozzi in [21] (see also references to previous partial solutions 
in jSniED])- 

In the flat situation, an amazing case was discovered by M. Avel- 
laneda and F.-H. Lin [Sj and later was also studied by J. Moser and 
M. Struwe jnHI- In these papers the authors dealt with polynomially 
growing solutions of second-order divergence form elliptic equations 

(1.3) Lu = - ^ {a'’^{x)ua;i)xj = 0 

with real coefficients that are periodic with respect to the lattice 7/^ in 
They obtained a comprehensive answer for (El (see also [221 E2] 
for related results and references). Using the formalism of homogeniza¬ 
tion theory, it was proved that the space of all solutions of the equation 
Lu = 0 of polynomial growth of order at most N has the same dimen¬ 
sion hn,N (see (US)) as the space of harmonic polynomials in of the 
same rate of growth. Moreover, any solution v of the equation Lv = 0 
in M"" of polynomial growth is representable as a finite sum of the form 


^We will also use the notation 

(1.1) 9".»^=("Y) 

for the dimension of the space of all polynomials of degree at most N inn variables. 
Notice that qn-i,N also coincides with the dimension of the space of all homogeneous 
polynomials of degree in n variables, so in particular, hn,N = <ln,N — <ln,N-2 = 
+ qn-l,N- 

^Without a condition on the curvature, the hyperbolic plane, where there is an in¬ 
finite dimensional space of bounded harmonic functions, provides a counterexample. 
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(1.4) v{x) = x^pj{x), 

where the functions Pj{x) are periodic with respect to the group of 
periods of the equation and x^ = 

One can say that there has been no complete understanding of this 
result concerning periodic equations. In particular, one can ask the 
following natural questions: 

(i) Is it important that the operator is of divergence form? 

(ii) Can the results be generalized to higher order equations? 

(iii) Is it possible to determine for a given periodic elliptic equation 
whether the Liouville theorem holds? 

(iv) How crucial is the usage of homogenization theory tools (which 
automatically restricts the class of equations)? 

(v) Same questions about elliptic systems. 

(vi) Can these results be generalized for covering spaces of compact 
manifolds? 

Some partial answers to these questions were obtained in |581 IKT] . In 
inn, the results for the divergence type operators (El in M"" were gen¬ 
eralized to the case of second-order periodic operators without lower 
order terms. At the same time, |2HI contains a necessary and sufficient 
condition for the validity of the Liouville theorem for a general periodic 
elliptic operator in as well as a description (in most cases, implicit) 
of the dimensions of the corresponding spaces of solutions. In partic¬ 
ular, an explicit formula was given for a general second-order periodic 
elliptic operator that admits a global positive solution. 

Simultaneously, Liouville type theorems for holomorphic functions 
on complex analytic manifolds were studied (see [ninniEniinniinniinni 
and references therein). For instance, one asks whether Liouville theo¬ 
rems for holomorphic functions on coverings of compact analytic man¬ 
ifolds (or more generally, on coverings of manifolds with the Liouville 
property) hold true. One should mention the results of 1221, where it 
was shown in particular that nilpotent coverings of compact complex 
analytic manifolds do have the Liouville property for hounded holomor¬ 
phic functions (i.e., the space of such functions is hnite-dimensional). 
It was not clear whether in these cases the spaces of holomorphic func¬ 
tions of a given polynomial growth are of hnite dimensions as well. 
The exception was the result of ra. where such a Liouville theorem 
was proven for abelian coverings of compact Kahler manifolds (see also 
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[2D]). This result also follows from EH, while this relation with Liou- 
ville theorems for harmonic functions disappears for non-Kahler case. 

One should also mention parallel studies concerning Liouville theo¬ 
rems for harmonic functions on graphs (e.g., j^IDH))- 

The goal of this paper is to provide results of Liouville type that clar¬ 
ify this issue for abelian covers of compact manifolds. The results apply 
to elliptic equations and systems (including overdetermined ones) on 
abelian coverings of compact Riemannian manifolds, as well as to holo- 
morphic functions on abelian coverings of compact complex manifolds, 
and to periodic equations on abelian coverings of combinatorial and 
quantum graphs. The crucial techniques used in the paper are differ¬ 
ent from the ones used in all the works cited above, with the exception 
of the authors’ paper |DH| • The ideas and the techniques come from the 
Floquet theory [snizzi and are related to some spectral notions com¬ 
mon in the solid state physics [7|. The reader can hnd all the necessary 
preliminary information in the next two sections. In comparison with 
[DHj, the present paper provides explicit formulas for the dimensions 
of the corresponding spaces, where (DHj, in general, contains only an 
implicit algorithm to calculate these numbers. Moreover, multiplicities 
at spectral edges, as well as overdetermined systems (in particular, d- 
operators and Liouville theorems for holomorphic functions) are now 
allowed. Furthermore, operators on graphs are also considered. The 
general theorems are applied to a variety of specihc examples of oper¬ 
ators. 

In order to outline the results of the paper, let us introduce some ob¬ 
jects. Consider a normal abelian covering^ of a compact d-dimensional 
Riemannian manifold M 

where G is the (abelian) deck group of the covering. Without loss of 
generality, one can assume that G = Z”. In fact, no harm will be done 
if the reader imagines for simplicity that X = M"', G = Z”, and M 
is the torus M^/Z” (albeit in general, the dimension d of M does not 
have to be equal to n). 

We will need to consider characters x of G, i.e., homomorphisms 
of G into the multiplicative group C* of nonzero complex numbers. 
Unitary characters map G into the group of complex numbers of 
absolute value 1. For any character y, a function f on X will be called 
X-automorphic if f{gx) = for any x E X, g E G. 

Let P be an elliptic G-periodic operator on X (in what follows, we 
use for shortness the word ‘periodic’ instead of ‘G-periodic’). For any 

^The word ‘covering’ in this paper always means ‘normal covering’. 
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character a, consider the space of x-aufoLnorphic functions on X. It 
can also be interpreted as the space of sections of a linear bundle over 
M determined by x- If is invariant with respect to P, so one can 
consider the restriction P(x) of P to this space"^. In the particular case 
of x{9) = 1) -P(l) is just the elliptic operator Pm on M, whose lifting to 
X is P. In ‘non-pathological’ cases, the spectra of all operators P(a) 
are discrete. The spectrum of P(a), as a multiple-valued function of 
the character y, is said to be the dispersion curve or dispersion relation. 

Let iV > 0. We say that the Liouville theorem of order N for the 
equation Pu = 0 holds true if the space Vn^P) of solutions of the 
equation Pu = 0 on X that satisfy |M(a;)| < C{l + p{x))^ for all a: G X 
is of hnite dimension. Here p{x) is the distance of a; G X from a hxed 
point Xo G X. 

We can now formulate a general (and somewhat vague at this point) 
statement that outlines our main results for the elliptic case contained 
in theorems in and M The results for the overdetermined, holomor- 
phic, and graph cases can be found in Sections El and [71 

Main Theorem. 

(1) If the Liouville theorem of an order N > 0 for the equation 
Pu = 0 holds true, then it holds for any order. 

(2) In order for the Liouville theorem to hold, it is necessary and 
sufficient that the number of unitary characters x for which the 
equation Pu = 0 has a nonzero x-outomorphic solution is finite. 

(3) If the Liouville theorem holds and the spectra of all the operators 
P{x) ore discrete, then under some genericity conditions on the 
operator P, the dimension of the space Vn^P) can be computed 
in terms of the dispersion relation for P. 

(4) If the Liouville theorem holds, then all solutions that belong to 
VAr(P) are linear combinations of Floquet solutions of order N 
(see the definitions in the text below). 

(5) Under the same conditions, one can describe a constant coeffi¬ 
cient (‘homogenized’) linear differential operator A(P) on M”, 
such that there is a one-to-one correspondence between the poly¬ 
nomial solutions of Av = 0 on M” and the polynomially growing 
solutions of Pu = 0 on X. 

We will see that this result in particular means that one should 
naturally expect the Liouville theorem to hold only when zero is at an 


^Exact definitions of function spaces and operators are provided in the next 
section. 
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edge of the spectrum of P. This is true, for instance, for the operators 
considered in jSl EH EH], when zero is the bottom of the spectrum. 

It is interesting to notice that the dimension of X does not have 
to be equal to n, so the operators P and A might act on manifolds 
of different dimensions. This happens since the Liouville property is 
of a ‘homogenized’ nature, i.e., it is something one sees by looking at 
the manifold X ‘from afar’. Thus, the local details of the manifold 
are essentially lost and one sees the euclidian space instead. In 
other words, the free rank of the deck group of the covering plays a 
more prominent role for Liouville theorems than the dimension of the 
manifold®. 

Similar results hold for elliptic systems, including overdetermined 
ones that are elliptic in the sense of being a part of an elliptic complex 
of operators. The reader can hnd basic notions and results concern¬ 
ing elliptic complexes in many books and articles (e.g., [TTl Vol. Ill, 
Section 19.4], [TH] Section 3.2.3], [HH Section IV.5]). For the partic¬ 
ular case of the Cauchy-Riemann d operator, one obtains a Liouville 
theorem for holomorphic functions on abelian coverings (Theorem EOj)- 
Analogs for operators on combinatorial and quantum graphs are also 
straightforward to obtain. 

The outline of the paper is as follows. The next section introduces 
necessary notations and preliminary results from the Floquet theory, in 
particular the dehnition and properties of the Floquet-Gelfand trans¬ 
form. The proofs are the same as for the case of periodic operators on 
and are hence mostly omitted (see, e.g., how they can be worked out 
in parallel to the flat case in |1H1 and also in 12111221 EHl). The crucial 
SectionEjis devoted to the detailed study of the so called Floquet-Bloch 
solutions. In Section E) we derive Liouville type theorems for elliptic 
systems. Section El provides some examples of applications to specihc 
periodic operators, and Section El treats overdetermined systems, in¬ 
cluding the case of analytic functions on complex manifolds. Graphs 
are briefly considered in Section[71 The last sections contain concluding 
remarks and acknowledgments. 

2. Notations and preliminary results 

Let us introduce hrst some standard notions of Floquet theory (see 
121 EH mi ), which we will adjust to the case of abelian covers (this 
does not require any change in the substance). 


®This resonates with M. Gromov’s notion of quasi-isometry, when the space of a 
covering might be indistinguishable from its deck group IMIISH] 
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Let X be a noncompact smooth Riemannian manifold of dimension 
d equipped with an isometric, properly discontinuous, and free action 
of a finitely generated abelian gronp G. The action of an element 
g ^ G on X E X will be denoted by gx. Consider the orbit space 
M = X/G, which due to our conditions is a Riemannian manifold of 
its own. We will assume that M is compact. Hence, we are dealing 
with an abelian covering of a compact manifold 

(2.1) TT : X ^ M{= X/G). 

Switching to a subcovering X ^ M (—;> M) with a compact M, one can 
eliminate the torsion part of G. In what follows, we could substitute 
M for M and hence reduce the group G to Z”' with some n G N. We 
will therefore assume from now on that G = 17^. This will not 
reduce the generality of the results. 

Let P be an elliptic operator of order m on X with smooth coef- 
hcients® that commutes with the action of G. Such an operator can 
be pushed down to an elliptic operator Pm on M (or conversely, P 
is the lifting of Pm to X). The ellipticity is understood in the sense 
of the nonvanishing of the principal symbol of the operator P on the 
cotangent bundle (with the zero section removed) T*X \ (X x {0}). 
The dual operator (the formal adjoint) P* has similar properties (in 
particular, P* is also G-periodic). Here the duality is provided by the 
bilinear rather than the sesquilinear^ form 

<gj >= j f{x)g{x) dx. 

X 

All the preparatory facts and main statements here hold for lin¬ 
ear periodic matrix operators that are either standard elliptic (some 
times called elliptic in the Petrovsky sense) or elliptic in the Douglis- 
Nirenberg sense (e.g., (23, im Vol. HI, Section 19.5], and jTHl Section 
3.1.2.1]). The only difference in the proofs between the scalar and 
system cases arises in the necessity of introducing appropriate spaces 
of vector-valued functions, exactly as it was done in m Section 3.4]. 
Doing this, however, would on one hand be very routine, and on the 
other hand would make reading the text more difficult. Bearing this in 
mind, we will provide detailed considerations for scalar linear elliptic 
operators only. They transfer with no effort to systems. 


®The smoothness condition can be significantly reduced (see the corresponding 
remarks in Section Ell. 

'^This is not essential, but simplifies somewhat the calculations. 
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For any quasimomentum A; G C" we denote by 7 ^ the character 
of G = Z” dehned as 7 ^( 5 ') = Here g = (( 71 , G and 

k-g = kigi+- ■ ■+kngn- We will also use the notation |( 7 1 = \gi\ + - ■ ■+\gn\ 
and for a multi-index j = (ji,..., jn) £ Z" we denote ( 7 -^ = Ylidi"- If 
k G the corresponding character is unitary. Due to the obvious 27r- 
periodicity of 7 ^ with respect to k, it is sufficient to restrict ourselves 
to the real vectors k in the Brillouin zone B = [—vr, vr]”, which is 
a fundamental domain of the reciprocal (dual) lattice G* = (27rZ)"'. 
Periodizing B (i.e., considering the torus = M"'/27rZ"'), one obtains 
the dual group T” to G = Z"'. 

For any k G C”, we define the subspace Ll{X) of LlciX) con¬ 
sisting of all functions /(x) that are y^-automorphic, i.e., such that 
f{gx) = 7 fc(( 7 )/(x) = e'^'®/(x) for a.e. x E X. Alternatively, this space 
can be defined as follows. We can identify 7 ^ with a one-dimensional 
representation of G and consider the one-dimensional flat vector bun¬ 
dle Ek over M associated with this representation. Then elements of 
L\{X) can be naturally identified with L^-sections of E^. A similar con¬ 
struction works also for other classes of functions (e.g., from Sobolev 
spaces). We will identify G-periodic (i.e., 70 -automorphic) functions on 
X with functions on M. Due to the periodicity of the operator P, it 
leaves the spaces Ll invariant, and so its restrictions to these subspaces 
define elliptic operators P{k) on the spaces of sections of the bundles 
Ek over M If P is selfadjoint, then P{k) is selfadjoint for any real 
quasimomentum k. 

It is natural that Fourier transform with respect to the periodicity 
group G reduces the space L^{X), as well as the original operator P 
on X, to the direct integral of operators P{k) on sections of Ek\ 


( 2 . 2 ) 


and 


© 

= j LHx)dk, 

B 


© 

(2.3) P = j P{k)dk. 

B 

The integral is understood with respect to the normalized Haar measure 
on the dual group T”, which boils down to the normalized Lebesgue 
measure dk on the Brillouin zone B. The isomorphism (in fact, an 

®In the case when X = K" with the natural Z” action, these operators can be 
identified with the “shifted” versions P{x,D + k) of the operator P acting on the 
torus T" = (see [HIEH)- 
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isometry) in (j2.2j) is provided by an analog of the Fourier transform 
(see jSH Section 2.2], [TTIIHSI), which we will call the Floquet-Gelfand 
transform U : 

(2.4) f{x) ^ Uf{k, a;) = ^ figxh-kig) k G C”. 

g&G 

This transform is the main tool in the Floquet theory for PDFs (e.g., 
isi [771 EH uni). It was introduced hrst in ina in order to obtain ex¬ 
pansions into Bloch generalized eigenfunctions for periodic selfadjoint 
elliptic operators. 

It is not hard to describe the image of a Sobolev space H^{X) under 
the Floquet-Gelfand transform. In order to do so, let us consider a 
quasimomentum k & and denote by the closed subspace of the 
space consisting of 7 fc-automorphic functions. It is clear that 

can be naturally equipped with the structure of a Hilbert space and 
that it can be identihed with the space H^{Ek) of if'^-sections of the 
bundle Ek over M. 

One can show® that 

(2.6) f := U Ht 

fceC" 

forms a holomorphic 27rZ"^-periodic Banach vector bundle. As any 
inhnite dimensional analytic Hilbert bundle over a Stein domain, it is 
trivializable [221 (see also the survey jHH] and theorems 1.3.2, 1.3.3, and 
1.5.23 in mi). 

We collect now several statements from [771 Theorem XHI.97], jHIl 
Theorem 2.2.2], and (2211221, recasted into the abelian covering form: 

Theorem 1. (1) For any nonnegative integer m, the operator 

U : ii™(X) ^ L2(T^,^™) 

is an isometric isomorphism, where L^(T”,T"*) is the space 
of square integrable sections over the torus (identified with the 
Brillouin zone F>) of the bundle equipped with the natural 
topology of a Hilbert space. 

(2) Let K (£ X be a domain in X such that U gK = X. Let also 

g£G 

the space 

sup||.).||„„|,71+|9|)"<«D Vivj 

®See Theorem 2.2.1 in ini for the case X = K”. The general case of abelian 
covers over compact manifolds is entirely parallel. 
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be equipped with the natural Frechet topology. Then 
U : C^{X) 

is a topological isomorphism, where is the space 

of sections of the bundle over the complex torus T”, 
equipped with the standard topology. 

(3) Let the elliptic operator P be of order m. Then under the trans¬ 
form U the operator 

P : C™(X) ^ C°(X) 

becomes the operator 






of multiplication by the holomorphic Fredholm morphism P{k) 
between the fiber bundles and E^. 


3. Floquet-Bloch solutions 

We now need to introduce and study our main notions: Bloch and 
Floquet solutions of periodic differential equations. 

Definition 2. Let k G C”. A q^-automorphic function u{x) on X is 
said to be a Bloch function with quasimomentum k. In other words, it 
is a function with the property u{gx) = 'yk{g)u{x) = e‘^'%(a;) for any 
X E X, g E G. To put it differently, u{x) is transformed according to 
an irreducible representation of the group G with the character 

A Bloch solution of an equation is a solution that is a Bloch function. 

Notice that every continuous Bloch function on X with a real quasi¬ 
momentum (i.e., transformed according to an irreducible unitary rep¬ 
resentation) is bounded. Any such Bloch function u that belongs to 
Li^^{X) is bounded in the following integral sense: for any compact 
K (£ X we have sup^g^ ||'w||L 2 (c,ic) < oo. 

In the case when X = R”, G = Z", and M = T", any Bloch function 
with quasimomentum k has the form 

u{x) = e'^'^p{x) 

with a Z"^-periodic function p{x). In fact, a similar (albeit less natural) 

representation holds for Bloch functions on any abelian cover X i—> M. 
Indeed, let K be any fundamental domain of X with respect to the 
action of G and / G (^“(X) be a nonnegative function which is strictly 
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positive on K. We define for any j = 1,... ,n 

hj{x) := fi9x)exp{-gj). 

9=(9i,---,5n)eG=Z" 

Then hj{x) clearly is a positive function satisfying hj{gx) = e^^hj{x) 
for any g = {gi,..., gn) G G. It is an analog of on Thus, one 
can define analogs of powers x^ = Xi ■ ■ ■x’’^ and of exponents e'^'^ as 
follows: for / = (/i, G Z" let 

n 

j=i 

and for any quasimomentum G C"' 

ek{x) := exp{iki\oghi{x) H-h iA:„ log 

Notice that ek{x) is a nonvanishing Bloch function on X with quasi- 
momentum k, which is positive for ik G R”. Thus, any Bloch function 
M on X with a quasimomentum k is given by 

u{x) = ek{x)p{x), 

where p{x) is G-periodic. 

The construction of Bloch functions can be described in a more in¬ 
variant way [2j. Consider a basis coj of the space of all closed differential 
1-forms on M (modulo the exact ones) such that their lifts Wj to X are 
exact. According to De Rham’s theorem, this basis is finite. One can 

X 

now achieve the same goals as before defining hj{x) = exp(f wj) for a 

O 

fixed reference point o E X. 

We can now define a more general class than Bloch functions. 


Definition 3. A function u{x) on X is said to be a Floquet function 
with quasimomentum fc G C” if it can be represented in the form 


/ 


\ 


(3.1) 


u{x) = ek{x) 



, 

\j\ <N 


where the functions pj are G-periodic. 

The number N in this representation will be called the order of the 
Floquet function. 

A Floquet solution of an equation is a solution that is a Floquet 
function. 
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This definition is modelled closely after the notion of Floquet solution 
that is common in M"" (e.g., ED), where the formula is the same, one 
just replaces the ‘powers’ [x]^ by the true powers xh However, unlike 
the notion of a Bloch function, it lacks invariance. To alleviate this, 
we briefly address now a different way to define Bloch and Floquet 
functions (solutions) on abelian coverings [fidj . 

For any g ^ G, let us define the first difference operator Ag acting 
on functions on X as follows: 

(3.2) Agu{x) = u{gx) — u{x). 

Clearly, u is a periodic function (i.e., a Bloch function with zero quasi¬ 
momentum) if an only if it is annihilated by Ag for any g E G. In fact, 
it is sufficient to check this property for any set {gj} of generators of G. 
One wonders whether one can check in a similar way whether a func¬ 
tion is a Bloch function with a nonzero quasimomentum and whether 
Floquet functions allow for similar tests. In order to get the answer, we 
need to introduce a twisted version of the first difference, that depends 
of the quasimomentum: 

(3.3) Ag.ku{x) = x-k{g)u{gx) - u{x) = e~'^'^u{gx) - u{x). 

We also need to introduce iterated finite differences of order N with 
quasimomentum k as follows: 

(3.4) Ag^^,,,^gj^-k = ^gi\k • . . . ■ Ag^-k, 

where gj G G (it will always be sufficient to use only elements (maybe 
repeated) of a fixed set of generators of G). 

We can now answer the question by proving the following 

Lemma 4. A function u{x) on X is a Floquet function of order N 
with quasimomentum k if and only if it is annihilated by any difference 
of order N + 1 with quasimomentum k: 

b ■ ■ ■ 1 9n+i E G, 

(choosing gj from a fixed set of generators is sufficient). In particular, 
it is a Floquet function of order N with quasimomentum 0 if and only 

(3.5) Ag^- ...■ = 0 V^fi,..., G G. 

Proof: Let us provide the proof for the case k = 0 first. As it has 
already been mentioned, the necessity and sufficiency of the condition 
(USD checks out easily for N = 0, where it boils down to AgU = 0 for 
all g E G, i.e., to the periodicity of u. Necessity for any N now follows 
easily by induction if one takes into account the representation (ED). 
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Indeed, one computes that (using the same standard basis {gj} of Z" 
as before) 

(3.6) + lower order terms. 

Here ‘lower order terms’ contain linear combinations of [a:]*s of strictly 
lower total degrees. Since the difference operators do not alter periodic 
functions, we obtain that any reduces Floquet functions of order 
A + 1 to the ones of order A, which concludes the induction step of 
the proof of necessity for fc = 0. 

Let us prove sufficiency which also follows by induction with respect 
to the order of the Floquet function. We have already checked it for 
A = 0. Assume that this has been proven for orders up to A. Suppose 
that u satishes A^^ -.. .-Ag^^^u = 0 for any ^fi,..., gN +2 ^ G- Take a set 
of generators {gj}'j=i in G. Then the functions fj := Ag.u satisfy the 
condition of the Lemma for the order A. According to the induction 
hypothesis, we conclude that 

/ \ 

I = (li, . . . ,ln) E Z” 

V \1\<N 

with periodic functions pij. 

We claim that for any Floquet function / of order A (with A; = 0) 
there exists a Floquet function v of order A + 1 such that Ag.v = f. 
Without loss of generality we assume that j = 1. Since the difference 
operator A^^ does not change periodic coefficients, it is sufficient to 
check the statement for / = [x\\ where |/| < A. Now induction with 
respect to li hnishes the job. Namely, (EH) provides a linear system of 
a triangular structure for recursively determining the coefficients of v 
such that Ag.v = /. 

Therefore, a Floquet function vi of order A + 1 exists such that 
AgjTi = AgjM. This means that the function U 2 = u — Vi is periodic 
with respect to the one-parameter subgroup generated by gi. Since 
the condition of the Lemma is still satished for M 2 , and since U 2 is gi- 
periodic, we can continue this process and hnd a Floquet function V 2 
of order A -|- 1 that is ^^i-periodic and such that uz = u — vi — V 2 is 
periodic with respect to both gi and g 2 . Continuing this process, we 
get the conclusion of the Lemma for fc = 0. 

Let us now prove the statement for an arbitrary quasimomentum 
k. According to the dehnition (EH) of Floquet functions, any Floquet 


fjA) = 
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function u of order N with quasimomentum k has the form u = ek{x)v, 
where u is a Floquet function of order N with quasimomentum k = 0. 
A straightforward calculation shows that if any two functions u and 
V satisfy u = ek{x)v, then one has Ag-^u = ^k{x)AgV. Thus, the 
statement of the lemma for an arbitrary quasimomentum k follows 
from the one for fc = 0. □ 

We have provided several different ways to interpret the notion of 
Floquet solutions: using explicit formulas analogous to the ones in R” 
(constructing analogs of coordinate functions and exponents by either 
explicit constructions, or by using some special differential forms on 
M), as well as in terms of some difference operators^®. Another way to 
think of Floquet functions is to imagine the hnite dimensional subspace 
E generated by the G-shifts of such a function as a Jordan block for 
the action of G [Hj. This suggests a relation to indecomposable (non¬ 
unitary) representations of G (e.g., [3 Ch. 6, Sect. 3]). However, we 
do not pursue this approach, due to the known difficulties of classifying 
such representations jSSj. 

Remark 5. Any continuous Floquet function u{x) of order N with a 
real quasimomentum satisfies the growth estimates 

|w(a:)| < G{1 + \g\)^ \/g & G and x G gK, 

where G depends on K and u. In general, one needs to replace this 
growth estimate by an integral one, as it was done before for Bloch 
functions. 

We now introduce a notion that plays a very important role in solid 
state physics, photonic crystal theory, as well as in the general theory 
of periodic PDFs HlEIlEaiZI]. It will be also crucial for formulation 
of our main results. 

Definition 6. The (complex) Fermi surface Fp of the operator P (at 
the zero energy level) consists of all quasimomenta k & such that 
the equation Pu = 0 on W has a nonzero Bloch solution with a quasi¬ 
momentum k . The real Fermi surface Fpr is Fp fl R"". 

Equivalently, k & Fp means the existence of a nonzero solution u of 
the equation P{k)u = 0. In the Euclidean case, such a u has a form 
u{x) = e^^'^p{x), where p{x) is a G-periodic function. Fermi surface 
plays in the periodic situation the role of the characteristic variety (the 
set of zeros of the symbol) of a constant coefficient differential operator. 

^^Difference operators approach has been successfully used in related studies 
of periodic equations and Liouville type problems for holomorphic functions in 

[HiniiMiini- 
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Introducing a spectral parameter A, one arrives at the notion of the 
Bloch variety. 

Definition 7. The (complex) Bloch variety Bp of the operator P con¬ 
sists of all pairs {k,X) G such that the equation Pu = Xu has 

a nonzero Bloch solution u with a quasimomentum k. The real Bloch 
variety Bp^^ is Bp fl 

The Bloch variety Bp can be treated as the graph of a (multivalued) 
function A(/c), which is called the dispersion relation. If the spectra 
of the operators P{k) on M are discrete, we can single out continuous 
branches Xj of this multivalued dispersion relation. They are called the 
band functions [ZZIEI] The Fermi surface is obviously the zero level 
set of the dispersion relation. 

In order to justify the notion of a band function, we need to guarantee 
the discreteness of the spectrum of the operators P{k) on M for all 
fc G C”. In other words, we need to exclude the pathological (but 
possible) situation of the spectrum of P covering the complex plane. 
One of the reasons for such strange spectral behavior is the Fredholm 
index being not equal to zero. However, even when the index is equal to 
zero, such pathology can occur, as the example of the operator e“^d/d(p 
on the unit circle shows. Self-adjointness of P is one of the conditions 
that would obviously guarantee discreteness. Another example is a 
second-order elliptic periodic operator in M” of the form 

n n 

(3.7) ^ “ X] aij{x)didj + ^ hi{x)di P c{x) 

i,j=l i=l 

with real and smooth coefficients. More sufficient conditions can be 
found for example in pp. 

Lemma 8. The Fermi and Bloch varieties are the sets of all zeros of 
entire functions of a finite order in and respectively. 

This is proven in m theorems 3.1.7 and 4.4.2] for the flat case. The 
case of a general abelian covering does not require any change in the 
proof. 

The lemma implies in particular, that the band functions X{k) are 
piecewise-analytic (e.g., when there is no level crossing, one has ana- 
lyticity due to the standard perturbation theory 01j). This statement 
was originally proven in jHS] for Schrodinger operators. 

Another useful property of the Bloch and Floquet varieties is the 
relation between the corresponding varieties of the operators P and 

p*. 
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Lemma 9. m Theorem 3.1.5] A quasimomentum k belongs to Fp* 
if and only if —k G Fp. Analogously, (fc, A) G Bp* if and only if 
(—fc, A) G Bp. In other words, the dispersion relations X{k) and X*{k) 
for the operators P and P* are related as follows: 


(3.8) 


y(k) = \{-k). 


We will need to see how the structure of the functions of Floquet type 
(see Dehnitionini) and in particular, of Floquet solutions of our periodic 
equation reacts to the Floquet-Gelfand transform. For instance, in the 
constant coefficient case, where the role of the Floquet solutions is 
played by the exponential polynomials 


u{x) = ^ 1 


lil<Y 


where pj G C, such functions are Fourier transformed into distributions 
supported at the point {—k). The next statement shows that under 
the Floquet-Gelfand transform, each Floquet type function (ED) cor¬ 
responds, in a similar way, to a (vector valued) distribution supported 
at the quasimomentum {—k). This, and some other properties of Flo¬ 
quet solutions that play a crucial role in establishing our Liouville type 
theorems are collected in the next lemma. 

Every Floquet type function u with a real quasimomentum is of 
polynomial growth, and thus determines a (continuous linear) func¬ 
tional on the previously defined space G°(X) (see Theorem P). If it 
satisfies the equation Pu = 0 for a periodic elliptic operator of order 
m, then as such a functional it is orthogonal to the range of the dual 
operator P* : G™'(X) —>• G°(X). According to Theorem [TJ after the 
Floquet-Gelfand transform any such functional becomes a functional 
on C°^ (T*^, £^) that is orthogonal to the range of the operator of mul¬ 
tiplication by the Fredholm morphism P*{k) : T"* ^ S^. The following 
auxiliary result (see IHH]) describes all such functionals^^. 

Lemma 10. (1) A continuous linear functional u on G°(X) is gen¬ 

erated by a Floquet type function with a quasimomentum ko if 
and only if after the Floquet-Gelfand transform it corresponds 
to a functional on C°^ (T’^,T°) which is a distribution cf that is 
supported at the point —k^, i.e., has the form 


(3.9) 



^^Although the discreteness of the spectrum of P{k) was assumed throughout 
the whole paper m, it is in fact not needed for (and was not used in the proof of) 
this lemma. 
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where qj G L‘^{M). The orders N of the Floquet function and 
of the corresponding distribution cf are the same. 

(2) Let Ok be the dimension of the kernel of the operator 

P{k) : ^ Ll 

Then the dimension of the space of Floquet solutions of the equa¬ 
tion Pu = 0 0 / order at most N with a quasimomentum k is 
finite and does not exceed akqn,N- 

The estimate on the dimension given in the Lemma ITUl is very crude 
and can often be improved. The next result (Theorem provides 
in some cases exact values of these dimensions. This theorem is the 
crucial part of the proof of (Liouville) Theorem In order to state 
and prove it, we need to introduce some notions. 

First of all, the analytic Hilbert bundle 8"^ is locally trivial for any m. 
Since the previous Lemma shows that our interest in Floquet solutions 
is local with respect to the quasimomentum fc, we can trivialize the 
bundles and hence assume that the analytic families of operators P{k) 
and P*{k) act in a hxed Hilbert space. At this moment, we will need the 
additional condition that the spectra of these operators are discrete (see 
the corresponding discussion earlier in the text). Assume now that zero 
is an eigenvalue of the adjoint operator P*(—fco) : ^-ko ^—ko ix) 

(since under the conditions we imposed on the operator, its Fredholm 
index is zero, this means that the operator P(fco) has eigenvalue zero 
as well). Suppose further that the algebraic multiplicity of the null 
eigenvalue is equal to r. Consider a closed curve T in C separating 
0 from the rest of the spectrum of P*{—ko) and the corresponding 
(analytically depending on fc in a neighborhood of ko) r-dimensional 
spectral projector n(A;) for P*{—k). Let be an orthonormal 

basis of the spectral subspace of P*(—fco) corresponding to the point 
0 G C (i.e., the range of n(fco)). Some of our main results will be 
expressed in terms of the following r x r matrix function 

(3.10) X{kY = {ejYP*{kMk)ei). 

This matrix function is analytic with respect to fc in a neighborhood 
of ko. 

Remark 11. We would like to mention that in what follows, the re¬ 
sults will be invariant with respect to a multiplication of the matrix A 
from either side by an invertible matrix-function analytic in a neigh¬ 
borhood of ko- This means, in particular, that if r = 1, then X{k) can 
equivalently be chosen to be equal to the analytic branch around ko of 
the eigenvalue of P*{—k) that vanishes at ko. 
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Moreover, instead of {Ii{k)ei} in (jd.lOjl one can nse any holoniorphic 
basis ej{k) in the range of n(A;) near k = k^. One can also nse instead 
of {ej, •) any family {fj{k)Yj^i of holomorphic with respect to k (in a 
neighborhood of /cq) fnnctionals that is complete on the range of n(A;). 

Consider the Taylor expansion of X{k) aronnd the point ko into ho- 
mogeneons matrix-valned polynomials: 

(3.11) \{k) = Y,Uk-h)- 

l>0 

In this paper we will be mostly concerned with the hrst nonzero term 
Xig of the expansion. 

Definition 12. Let Q be a homogeneons polynomial in n variables with 
matrix coefficients of dimension r x r, and let Q{D) be the differential 
matrix operator with the symbol Q. A C^-valned polynomial p{x) in R” 
is called Q-harmonic if it satisfies the system of differential eqnations 
Q{D)p = 0. 

Let V denote the vector space of all C^-valned polynomials in n vari¬ 
ables, and let P; be the snbspace of all snch homogeneons polynomials 
of degree 1. So, 

OO 

(3.12) ^ = 0P/ 

1=0 

and 

N 

(3.13) P^;=0P, 

1=0 

is the snbspace of all snch vector valned polynomials of degree at most 
N. If Q{k) is a homogeneons polynomial of degree s with valnes in 
r X r matrices, then the matrix differential operator Q{D) maps Pi+s 
to P;. If the determinant det Q is not identically eqnal to zero, then 
this mapping is snrjective for any I (this will follow from the proof of 
the theorem below). Hence, the mapping Q{D) : V —>■ V has a (non- 
nniqnely dehned) linear right inverse R that preserves the homogeneity 
of polynomials. 

Theorem 13. Assume that zero is an isolated eigenvalue of algebraic 
multiplicity r of the operator P*{—kf) : iL!!\,^(X) —>• L^^^(X),. Let 
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\{k) be defined in a neighborhood of ko as in . Let also Az(, be 

the first nonzero term of the Taylor expansion \d.ll\) . Then 

(1) For any > 0, the dimension of the space of Floquet solutions 
of the equation Pu = t) in X with quasimomentum k^ and of 
order at most N is finite and does not exceed rqn,N- 

(2) If det A/q is not identically equal to zero (for instance, this is 
the case when the eigenvalue is simple, i.e., r = 1), then for 
any N > 0 the dimension of the space of Floquet solutions of 
the equation Pu = 0 in X of order at most N and with quasi¬ 
momentum ko is equal to 


(3.14) 


n + iVA / n-\- N — Iq 
N J ^ N-lo 


(3) This dimension coincides with the dimension of the space of all 
XiQ-harmonic polynomials of degree of at most N with values in 
C^. Moreover, given a linear right inverse R of the mapping 
XifiD) : V ^ V that preserves homogeneity, one can construct 
an explicit isomorphism between the corresponding spaces. 


Proof: As before, we will assume that the bundles are analytically 
trivialized around the point fco and hence all operators act between 
hxed spaces that we will denote and H^. In order to simplify 
notations, we assume that fco = 0 (this can always be achieved by a 
change of variables). Let us denote by N{k) the range of the projector 
n(A;) and choose a closed complementary subspace M to A^(0) in H"^. 
The subspace M stays complementary to N{k) in a neighborhood of 0 
and so 


= M®N{k). 

Thus, P*{—k) has zero kernel on M for all A; in a neighborhood of 0. 
This implies that the range of P*{—k) on M forms an analytic Banach 
vector bundle R{k) in a neighborhood of 0 (e.g.. Theorem 1.6.13 of |51jj. 
Representing now the operator P*{—k) in the block form according to 
the decompositions 

H^ = M® N{k) 

and 

= R{k)®N{k), 


we get 


P*{-k) 


B{k) 0 A 
0 X{-k) ) ’ 


^^Any analytic function in a neighborhood of ko that differs from X{k) by a left 
and right multiplication by analytic invertible matrix functions will produce the 
same results in what follows. 
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where B{k) is an analytic invertible operator-function and the matrix 
analytic function X{—k) differs from X{—k) only by multiplying by an 
invertible analytic matrix function and thus for our purposes can be 
replaced by the latter one. 

Let us now have a functional 0 on (T”, So) supported at 0, such 
that it is orthogonal to the range of the operator of multiplication by 
P*{—k). Then it must be equal to zero on all sections of the bundle 
R{k) (since they are all in the range, due to invertibility of B{k)). 
This means that the restriction of such functionals to the sections of 
the hnite-dimensional bundle N{k) is an one-to-one mapping. This 
reduces the problem to the following: hnd the dimension of the space 
of all distributions of order N supported at the origin such that they 
are orthogonal to the sub-module generated by the matrix X{—k) in 
the module of germs of analytic vector valued functions. One can 
change variables to eliminate the minus sign in front of k. Due to the 
hniteness of the order of the distribution, the problem further reduces 
to the following: hnd the dimension of the cokernel of the mapping 


Aat : Vn —^ Vn ■ 


Here An{p) for p &Vn is the Taylor matrix-valued polynomial of order 
TV at 0 of the product \{k)p{k). Let us write the block matrix Aj^ of 

N 

the operator Ajv that corresponds to the decomposition Vn = 0 Pp 

z=o 

Then Ajj = 0 for i — j < /q. For i — j > lo the entry A^ is the operator 
of multiplication by \i-j acting from Pj into P*. Since det(AzQ) is not 
identically equal to zero, it follows that for i — j = Iq the operator Aij 
of multiplication by Xi^ has zero kernel. 

In order to prove the theorem, we need to hnd the dimension of the 
cokernel of An, as well as to obtain the cokernel’s description. 

The hrst statement of the theorem is now obvious, since the dimen¬ 
sion of the cokernel of An cannot exceed the dimension of the ambient 
space Vn, which is equal to rqn^N- 

Let us now approach the second statement. Since Aat is a square 
matrix, dim Coker Aat = dim KerAAr. The latter dimension, however, 
is easy to hnd, due to the triangular structure of the equation AatP = 0, 
if it is written in the block matrix form according to the decomposition 
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(jini): 



/ 0 ... 0 A/q A/q+i ... Xn \ 


/ Pn \ 


0 ... 0 0 Aip ... Atv-i 


Pn-1 

(3.15) 

0 

<< 

•• 0 

•• 0 


PN-Io 


0 ... 0 


Pn-Io-1 


iyO ... 0 J 


\ Po J 


Here p = Po + ■ ■ ■ + Pn is the expansion of p G Vn into homogeneous 
terms. Since detAig 7 ^ 0, one concludes immediately from (Id.1511 that 
Po = ... = Pn-Io = 0, while other components are arbitrary. This gives 
the dimension of the kernel (and hence of the cokernel) of A at as 

r/n + iV\_/n + iV-/o Y 

[[ N J [ N-lo JJ’ 

which proves the second statement of the theorem. 

In order to prove the last assertion of the theorem, we need describe 
the elements of the cokernel of Ajv. Hence, we need to find the kernel 

N 

of the adjoint matrix AJ^. The adjoint matrix acts in the space 0P^, 

1=0 

where can be naturally identihed with the space of linear combina¬ 
tions with coefficients in of the derivatives of order I of the Dirac’s 
delta-function at the origin. Here we have A*^ = 0 for j — i < Iq, and for 
j — i > Iq the entry A*j is the dual to the operator of multiplication by 
Xj-i acting from Pj into Pj. In particular, since for j — i = Iq the latter 
operator is injective, we conclude that the operators A*j are surjective. 
This enables one to describe the structure of the kernel of A)^ (and 
hence of cokernel of Atv). Namely, let 

N 

V’ = ii’o, ■ ■ ■, i>N) G ^ P* 

1=0 

be such that A^i/) = 0. Due to the triangular structure of A^, we can 
solve this system: 

i = 0,...,N -k. 

j>i+lo 

Taking the Fourier transform, we can rewrite this system in the form 
Xj-i{D)'il:j = 0 i = 0,...,N -Iq, 

j>i+lo 
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where ip denotes the Fourier transform of ip. Therefore, ipj is a homo¬ 
geneous polynomial of degree j in M”. For i = N — Iq we have 

\iD)ipN = 0 . 

This equality means that ipN can be chosen as an arbitrary A^p-harmonic 
homogeneous polynomial of order N (with values in C^). Moving to 
the previous equation, we analogously obtain 

+ ^lo+l{D)lpN = 0 , 

or 

Mo{D)lpN-l = —XlQ+l{D)lpN. 

The right hand side is already determined, and the nonhomogeneous 
equation, as we concluded before, always has a solution, for instance 

-R (^Xig+i{D)ipN^ ■ 

This means that 

'0Y-1 + R ^Aip+i(D)'0vj 

is a A/p-harmonic homogeneous polynomial of order N — 1. We see that 

the solution ipN-i exists and is determined up to an addition of any 
homogeneous A^p-harmonic polynomial of degree — 1. Continuing 
this process until we reach ipQ, we conclude that the mapping 

Ip = {ipo, . . . , iPn) ^ (p= (00, •••,0a), 

where 

<pj ipj -|- -R ^ ^ Xi—j^if^(^D')ipi 
i>j 

establishes an isomorphism between the cokernel of the mapping 
and the space of A^p-harmonic polynomials of degree at most N. This 
proves the theorem. □ 

In the simplest possible cases, the theorem immediately implies the 
following: 

Corollary 14. Under the hypotheses of TheoremM,‘t\ one has: 

(1) If ko is a noncritical point of a single band function X{k) such 
that A(0) = 0, then the dimension of the space of Floquet so¬ 
lutions of the equation Pu = 0 on X of order at most N with 
a quasimomentum ko is equal to qn-i,N, be., it is equal to the 
dimension of the space of all polynomials of degree at most N 
in 
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(2) If the Taylor expansion at a point ko of a single band function 
A(/c) such that A(0) = 0 starts with a nonzero quadratic formf^, 
then the dimension of the space of Floquet solutions of the equa¬ 
tion Pu = t) on X of order at most N with quasimomentum 
is equal to hn,N, it is equal to the dimension of the space 
of harmonic (in the standard sense) polynomials of degree at 
most N in M”. In particular, this condition is satisfied at non- 
degenerate extrema of dispersion curves, i.e., the condition is 
satisfied at nondeqenerate spectral edges. 

In both cases an isomorphism can be provided explicitly as in the pre¬ 
vious theorem. 


4. Liouville type theorems for elliptic periodic systems 


In this section we shall use the results of the previous sections to 
establish Liouville theorems for periodic equations. We will consider 
at the moment an arbitrary linear (square) matrix elliptic operator 
P on the space X of the abelian covering X ^ M with smooth G- 
periodic coefficients that satishes the assumptions made in Section 
As before, without loss of generality we can limit the consideration 
to the case G = IP. Any of the standard meanings of ellipticity of 
a system would do, e.g., ellipticity in Petrovsky or Douglis-Nirenberg 
sense [2H1- 


Definition 15. Let K (£ X he a domain in X such that 
For iV = 0,1,..., dehne the space 


u oa: 

g£G 


X. 


V n{P) •= 3 I Pu = 0 in X, and 


sup 

geGL 


l“llL2(gA)(l + l5'l) 


-TV 


< OO 


We say that the Liouville theorem of order N > 0 holds true for the 
operator P if dim VAr(P) < oo. 

We say that the Liouville theorem holds true for the operator P if it 
holds for any order N > 0. 


Abusing notations, we will call solutions in Vo(P) bounded solutions. 
Since obviously Vn^^P) C Vn^^P) for Ni < N 2 , Liouville theorem of 
higher order implies the lower order ones. It is not clear a priori that 
the converse holds. The next results show that this is in fact true, at 

^^This form can be degenerate, as opposed to our assumption in isa Corollary 
11 ]. 

^“^The results hold with the same proofs for periodic elliptic equations in sections 
of periodic vector bundles on X. 
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least in our situation of abelian coverings. This observation apparently 
has failed to be made in previous studies, which sometimes lead to 
investigations of some individual Liouville theorems, e.g., for iV = 0 
without noticing their simultaneous validity for all N. 

We will start with an auxiliary statement, which is an analog of the 
classical theorem on the structure of distributions supported at a single 
point. It is a generalization of Lemma 25 in |2H]; where Fredholm rather 
than semi-Fredholm property is assumed. For the results of this section 
Lemma 25 in would be sufficient, while we need the full strength 
of the lemma below to treat overdetermined systems and holomorphic 
functions in Section El 

Lemma 16. Let T he a -manifold and P : T ^ L{Hi, H 2 ) he a C°^- 
function with values in the space L{Hi, H 2 ) of hounded linear operators 
between the Hilhert spaces Hi and H 2 - Assume that for each k eT the 
operator P{k) is right semi-Fredholm (e.g., be., it has a closed 

range and a finite dimensional cokernel. Then 

(1) If P{k) is surjective for all points k in T, then the multiplication 
operator 

(4.1) C^{T, Hi) C^{T, H 2 ) 
is surjective. 

(2) For any fixed ko E T the dimension of the space of functionals 
of the form 

(4.2) 

that are orthogonal to the range of the multiplication operator 
o is finite. Here Qj are continuous linear functionals on H 2 , 
the pairing < gj^f) > denotes the duality between ^.nd H 2 , 
Dj^k are linear differential operators with respect to k on T, and 
N> 0. 

(3) If P{k) is surjective for all points k except of a finite subset 
F C T, then any continuous linear functional g on the space of 
smooth vector functions C°°{T, H 2 ) that annihilates the range 
of the multiplication operator 




J<N 




LIOUVILLE THEOREMS ON ABELIAN COVERINGS 


25 


has the form 

(4.4) <g^(j)>=Y^ 9j,(j)>) , 

ki&F \j<N J 

in the notations of the previous statement of the lemma. 

Proof: Let us establish the validity of the hrst statement of the lemma. 
Due to the existence of partitions of unity, the statement is local. Lo¬ 
cally, following for instance the proof of Theorem 2.7 in [SH], one can 
construct a smooth one-sided (right) inverse Q{k) to the operator func¬ 
tion P{k). Now multiplication by Q{k) provides a right inverse to (j4.1 jl . 
which proves the surjectivity. 

To prove the second statement, let us consider a closed subspace 
M G Hi complementary to the kernel of P{kQ). Then the operator 
P{kii) : M ^ H 2 is injective and Fredholm. This injectivity property 
is preserved in a neighborhood U of ko- In particular, the subspace 
M{k) = P{k){M) C H2 of hnite codimension forms a smooth subbun¬ 
dle Ai in U X H 2 U. Now, any smooth 772-valued function f{k) such 
that f{k) G M{k), belongs to the range of the operator (14.111 . Hence, 
functionals orthogonal to the range can be pushed down to the smooth 

sections of the hnite dimensional bundle C = U H 2 /Mik) over U. It is 

keu 

clear that the functional one gets on this bundle preserves the structure 
(03). Such functionals on a hnite dimensional bundle, however, form 
a hnite dimensional space (for any hxed N). 

The third statement can be proven analogously to the similar state¬ 
ment in 011 Corollary 1.7.2]. Namely, under the conditions of the 
statement, and taking into account the hrst claim of the lemma, any 
functional annihilating the range of the operator of multiplication by 
P{k) must be supported at the hnite set F over which P{k) is not 
surjective. We can reduce the consideration to a neighborhood U oi a. 
single point fco ^ F- Now, the proof of the second statement reduces 
the functional to one supported at the point fco and dehned on smooth 
sections of a hnite-dimensional bundle C over U. Thus, the standard 
representation of distributions supported at a point implies (03). □ 

The next theorem shows that the existence of a polynomially growing 
solution implies the existence of a nonzero bounded Bloch solution (i.e., 
a solution automorphic with respect to a unitary character of G). 

Theorem 17. The equation Pu = 0 has a nonzero polynomially grow¬ 
ing solution if and only if it has a nonzero Bloch solution with a real 
quasimomentum (such a solution is automatically bounded), i.e., if and 
only if the real Fermi surface Fp^R = Fp DMT is not empty. 
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Proof; Assume that Fp^r = 0. Then P{k) is surjective for all k G T"". 
Indeed, if this were not the case, we could hnd a nonzero functional of 
the type (Hi) with = 0 and some k. According to the hrst state¬ 
ment of Lemma cni this would mean the existence of a Bloch solution. 
Now, the hrst statement of Lemma uni guarantees the surjectivity of 
the mapping 




P(fc) 


^oo(-j,n,^0) 


and hence the absence of any nontrivial functionals on S^) that 

annihilate the image of this mapping. Since under the Floquet-Gelfand 
transform U, any polynomially growing solution u{x) of Pu = 0 is 
mapped to such a functional, we conclude that u = 0. □ 

The next result provides necessary and sufficient conditions under 
which Liouville theorems hold for equations of the type we consider. 
It also establishes that the validity of Liouville type theorems does not 
depend on the order of polynomial growth. 


Theorem 18. Under the conditions we have imposed on the covering 
X and the operator P, the following statements are eguivalent: 

(1) The number of points in the real Fermi surface Pp^R is finite 
(i.e., Bloch (or automorphic) solutions exist for only finitely 
many unitary characters jk)- 

(2) There exists N > 0 such that the Liouville theorem of order N 
holds true. 

(3) The Liouville theorem holds (i.e., it holds for any order N). 

Proof: (2) ^ (1). Any Bloch solution with a real quasimomentum k 
(i.e., corresponding to a unitary character) is bounded and hence be¬ 
longs to the space VAr(P) for any N. Since such solutions with different 
characters are linearly independent, the validity of the Liouville theo¬ 
rem for some value of N implies that the number of the corresponding 
characters is hnite. Since the characters of all Bloch solutions with real 
quasimomenta constitute the real Fermi variety, Pp^R is hnite. 

The implication (3) ^ (2) is obvious. 

Let us now prove (1) ^ (3). Let u be a nonzero polynomially grow¬ 
ing solution. It can be interpreted as a continuous functional on C° 
annihilating the range of the dual operator P* : C"^ —C°. After the 
Floquet-Gelfand transform U, one obtains a functional on 
orthogonal to the range of the operator 




P{k) 
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By our assumption, is finite, thus the second and the third state¬ 
ments of Lemma nni and Theorem usi hnish the proof of statement 
(3). □ 

While this theorem establishes conditions under which Liouville the¬ 
orem holds, it does not tell much about the dimensions of the spaces 
Vn{P) of polynomially growing solutions, besides those being hnite. It 
also does not address the structure of these solutions. The next result 
provides some estimates and even explicit formulas for the dimensions, 
as well a representation for the solutions. 

Theorem 19. Suppose that the Liouville theorem holds for an elliptic 
operator P and let dN ■= dimVAr(P). Then the following statements 
hold: 

(1) Each solution u G VAr(P) can be represented as a finite sum of 
Floquet solutions: 


(4.5) u{x)= 

3 

where each Ukj is a Floquet solution with a quasimomentum k, 
and Pp^K = Pp n M”". 

(2) For all N eN, we have 


dN < doQn^N < OO , 


where qn,N is the dimension of the space of all polynomials of 
degree at most N in n variables. 

(3) Assume that the spectra of the operators P{k) are discrete for all 
k and that for each real quasimomentum k G Pp^R the conditions 
of Theorem^^ are satisfied. Then for each N > 0 the dimension 
dN of the space VAr(P) is equal to 


(4.6) 


51 


n + N 
N 


n + N - lo{k) \ 
N-lo{k) ) 


Here and lo{k) are respectively the multiplicity of the zero 
eigenvalue and the order of the first nonzero Taylor term of 
the dispersion relation at the point k G Pp^R (see these no¬ 
tions explained before Theorem f73|) . The terms in this sum are 
the dimensions of the spaces of -harmonic polynomials, 

and polynomially growing solutions can be described in terms of 
these polynomials analogously to Theorem.W^ 


Proof: All these statements follow immediately from Lemma The¬ 
orem [13 Lemma IM and Theorem ITHl □ 
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5. Examples oe Liouville theorems for specieic operators 

Let us recall that the L^-spectrum of the operator P is the union 
over fc G B of the spectra of P{k) (e.g., (HU Theorem 4.5.1] and [771 
Theorem XIII.85]). In other words, the spectrum of P coincides with 
the range of the dispersion relation over the Brillouin zone B. We 
have also discussed that the real Fermi surface for P is just the zero 
level set for the dispersion relation over the Brillouin zone. Theorem 
UHl of the previous section shows that the Liouville theorem holds if 
and only if this Fermi surface is finite^^. One expects this to happen 
normally at ‘extrema’ of the dispersion relation (albeit this is neither 
necessary, nor sufficient). In other words, speaking for instance about 
the selfadjoint case, one should expect Liouville theorem to hold mainly 
when zero is at the edge of a spectral gap, although it is possible in 
principle to have interior points of the spectrum where such a thing 
could occur as welh®. One notices that the cases considered in (HI 
irm ins] all correspond to zero being at the bottom of the spectrum 
(in the second-order nonselfadjoint case we mean by the bottom of 
the spectrum the generalized principal eigenvalue, see niHHiinii and 
Section 5.4 below). This explains why the homogenization techniques 
employed in these papers could be successful. Indeed, homogenization 
works exactly at the bottom of the spectrum. The results of this work 
show that one should also consider internal spectral edges, where the 
standard homogenization does not apply. 

So, let us now try to look at some examples where one can apply 
the results of the previous sections. Theorem ^1 is applicable to any 
elliptic periodic equation or system of equations on any abelian covering 
of a compact manifold. According to this theorem, one only needs to 
establish the finiteness of the real Fermi surface. On the other hand. 
Theorem m which provides formulas for the dimensions of the spaces 
of polynomially growing solutions, as well as some representations of 
such solutions, is more demanding. It requires that one guarantees 
discreteness of the spectra of the ‘cell’ operators P{k), and most of all, 
requires understanding of the analytic structure of the dispersion curve 
near the Fermi surface. Concerning this structure, it is expected that 
the following is true: 

Conjecture 20. Let P be a ‘generic’ selfadjoint second-order elliptic 
operator with periodic coefficients, and let (A_, A+) be a nontrivial gap 
in its spectrum. Then each of the gap’s endpoints is a unigue (modulo 

particular, outside the spectrum the Liouville theorem holds vacuously, 
according to Theorem d and Theorem 4.5.1 in m 

^®Two spectral zones with touching edges could provide such an example. 
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the dual lattice) and nondegenerate extremum of a single band function 

X,{k). 

This conjecture is crucial in many problems of mathematical physics, 
spectral theory, and homogenization (see its further discussion in Sec¬ 
tion ED and is widely believed to hold. Unfortunately, the only known 
theorem of this kind is the recent result of [IHl; which states that gener- 
ically a gap edge is an extremum of a single band function. A similar 
conjecture probably holds for equations on abelian coverings of com¬ 
pact manifolds. Theorem ITHl shows that validity of Coni ectureOHl would 
imply the following statement: 

Theorem 21. If Conjecture \2(A holds true, then ‘generically’, at the 
spectral edges of a selfadjoint elliptic operator of second-order with pe¬ 
riodic coefficients on the Liouville theorem holds, and for each 
> 0 the dimension of the space Vis egual to h^^N, the dimen¬ 
sion of the space of all harmonic polynomials of order at most N in n 
variables. 

In the similar situation on abelian coverings of compact manifolds, 
the rank of the deck group would appear in the answer rather than the 
dimension of the manifold (see Section ED- 

At the bottom of the spectrum, however, much more is known na 
[IHl[IHlEIlll3IZ3Eni- Even this limited (in terms of spectral location) 
information together with theorems^! and ^1 provides one with many 
specihc examples that go far beyond the equations considered in jaEB 

IB^] . 

The hrst trivial remark is that if zero is outside the spectrum of the 
operator P, then according to Theorem ITHI the Liouville property holds 
vacuously. Indeed, in this case the real Fermi surface is empty, and 
hence the equation Pu = 0 has no polynomially growing solutions. Let 
us look now at some less trivial examples. 

5.1. Schrodinger operators. Let X —> M be, as before, a non¬ 
compact abelian covering of a d-dimensional compact manifold and 
H = —A -|- V{x) be a Schrodinger operator on X with a periodic 
real valued potential V G r > d. Then the result of jlB] for 

A = and of jlH| in the general case states that the lowest band 
function Ai(fc) has a unique nondegenerate minimum Aq at k = 0. All 
other band functions are strictly greater than Aq. In particular, the 
bottom of the spectrum of the operator H is at Aq. Let us assume 
that Aq = 0 (or replace the operator hy H = —A + V{x) — Aq). Then 
theorems M and M become applicable, since the real Fermi surface 
consists of a single point k = 0 and the band function has a simple 
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nondegenerate minimum at this point (i.e., r = 1 and /q = 2 in the 
notations of Theorem uni)- Thus, one obtains 

Theorem 22. (1) The Liouville theorem holds for the operator H— 

Aq. 

(2) The dimension of the space VAr(iJ — Aq) equals hn,N (where n 
is the free rank of G). 

(3) Every solution u G Yn{,H — Aq) of the equation Hu — Aqu = 0 
is a Floquet solution of order N with quasimomentum k = 0, 
i.e., it can he represented as 

u{x)= [xYPjix) 

\j\<N 

with G-periodic functions pj. 

5.2. Magnetic Schrodinger operators. Although the first result of 
this subsection holds also in the general situation of abelian coverings, 
we will present it for simplicity in the flat case. 

Consider the following selfadjoint magnetic Schrodinger operator on 
M" 

(5.1) H =[iV + A{x)f + V{x) 

with periodic electric and magnetic potentials V and A, respectively^^. 
The introduction of a periodic magnetic potential into the operator is 
known to change properties of the Fermi surface significantly (e.g., m 
EO]). However, small magnetic potentials do not destroy the properties 
of our current interest. Indeed, let V be as in the previous statement, 
then there exists £ > 0 (depending on V) such that for any periodic 
real valued magnetic potential A such that 

I I A| I < £ 
j A{x) da; = 0 

the following holds true: The lowest band function Xi{k) of H attains 
a unique nondegenerate minimum A^g at a point ko (albeit, not neces¬ 
sarily at fc = 0). All other band functions are strictly greater than A^^,. 
Indeed, when both V and A are sufficiently small, this is proven in EH, 

^^One can see that the considerations and the result of this subsection immedi¬ 
ately extend to more general abelian coverings. We avoid doing so, in order not to 
introduce any new notions required for defining the magnetic operators in such a 
setting. 


and 

(5.2) 
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It is not hard, though, to allow for arbitrary electric and small mag¬ 
netic potential. Indeed, the case when A = 0 is covered by the result of 
(see the previous subsection). Now, the statement of Lemma|Hl(see 
also m Theorem 4.4.2]) can be easily extended without any change 
in the proof to include analyticity with respect to the potentials (e.g., 
[!TT]b Namely, there exists an entire function f{k,X,A,V) of all its 
arguments such that f{k,\,A,V) =0 is equivalent to 

(/c. A) G i?(iv+A)2+V , 

where as before is the Bloch variety of the operator H. This, 
together with the just mentioned result of na for A = 0, imply that 
the lowest band function Xi{k) of H attains a unique nondegenerate 
minimum at a point fco for sufficiently small magnetic potentials'^. 
Now one uses theorems UHl and m again to obtain 

Theorem 23. Assume that V, A, and k^ satisfy the above assump¬ 
tions, then 

(1) The Liouville theorem holds for u G Vn{H — A^p). 

(2) Any solution u G Vjv(i/ — A^p) is representable in the Floquet 
form 

v{x) = ^ x^pj{x) 

\3\<N 

with periodic functions pj {x). 

(3) The dimension of the space Vn^H — A^p) is equal to hn,N- 

Note that the normalization (Q can be always achieved by a gauge 
transformation, which does not affect the spectrum and the Liouville 
property. 

Besides magnetic Schrodinger operators with small magnetic poten¬ 
tials, another very special subclass of periodic magnetic Schrodinger 
operators is formed by the so called Pauli operators. Consider the 
following Pauli operators in 

(5.3) P± = (iV + A)2 ± (a,^A2 - a,,Ai), 

where A[x) = {Ai[x), A 2 {x)) is periodic. The structure of the dis¬ 
persion curves at the bottom of the spectrum of such an operator was 
studied in ra- It was shown that the dispersion relation for P± attains 
at the point /c = 0 its single nondegenerate minimum with the value 
Aq = 0. This implies that the result for the Pauli operators holds ex¬ 
actly like for the Schrodinger operator in the previous subsection, and 
in fact in a more precise form, since the minimal value Aq of the band 

^®The degree of smallness of A depends on the electric potential V. 
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function is known. Namely, every solution u G V 7 v(^^±) of the equation 
Pj-U = 0 on is representable in the form (Q. The dimension of the 
space V 7 v(T±) is equal to h 2 ^N = 2iV + 1, which is the dimension of the 
space of all harmonic polynomials of order at most N in two variables. 

The cases above of small magnetic potentials and of Pauli operators 
form a rather special subclass of periodic magnetic Schrodinger opera¬ 
tors, in the sense that one still hnds a single nondegenerate minimum 
of band functions at the bottom of the spectrum. However, this is 
not true anymore for the whole class of periodic magnetic Schrodinger 
operators, which should influence signihcantly our Liouville theorems. 
Although there is not much known here yet, the results of ISDI pro¬ 
vide the hrst glimpse at the possibilities. In that paper the following 
operator in is considered: 

Mt = {iV + tA{x)Y - 

where t G M and the magnetic potential A{x) is (0, a(xi)) with the 
1-periodic function a{xi) that is equal to sign(xi) on [—0.5, 0.5]. The 
main result of is that for \t\ < 2\/3 the bottom of the spectrum 
of Mt is 0 and it is attained at the quasimomentum fc = 0, where 
the band function has a simple nondegenerate minimum. This implies 
immediately the same Liouville theorem as for the Pauli operators. 
The situation changes for |f| = 2a/3, when this minimum (which is 
still 0 and is attained at fc = 0) becomes degenerate. Namely, it is 
shown in |HD] that the lowest eigenvalue Xi{k) of the operators Mtik) 
for |t| = 2^/3 has the following Taylor expansion at k = 0: 

This shows the degeneration of the Hessian. On the other hand, since 
the hrst nonzero homogeneous term of the expansion is of order two. 
Theorem cni still says that the Liouville theorem sounds the same as 
before (i.e., dim V 7 v(M^ 2 x/ 3 ) = 2iV -|- 1). This changes, though, for a 
little bit larger values of |f|, when the bottom of the spectrum shifts 
to the negative half-axis and two points of extremum appear with a 
nonvanishing Hessian at these quasimomenta. This implies that the 
Liouville theorem will apply after shifting to the bottom of the spec¬ 
trum (i.e., after adding an appropriate positive scalar term to Mt), but 
the dimensions of the spaces Vat of polynomially growing solutions will 
increase above h 2 ,N- 

Theorem 24. (1) The Liouville theorem holds for the operators P± 

and Mt for |t| < 2^/3. 
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(2) The dimensions of the corresponding spaces Vn are equal to 
h2,N — 2iV + 1 . 

(3) Every solution u eVn can be represented as 

'^(^) = x^pj{x) 

\j\<N 

with G-periodic functions pj. 

(4) For a sufficiently small £ > 0 and the values |t| G (2 a/3, 2a/3 + 
e), the Liouville theorem holds at the bottom of the spectrum of 
the operator Mt, while the dimension of the space Y]si{Mt) is 
equal to 2h2,N- 

It is also conjectured in IHHI that for some values of t, the operator Mt 
with the magnetic potential equal to A{x) = (sign(a; 2 ), sign(a;i)) on the 
cube [—0.5, 0.5]^ and Z^-periodic, will exhibit a complete degeneration 
of the bottom of the spectrum in the sense of complete vanishing of the 
Hessian. If this happens to be true, then according to Theorem it 
would lead to an increase in the dimensions of the spaces Vat. 

5.3. Operators admitting regular factorization. A very general 
and wide class of selfadjoint second-order periodic operators H on M”, 
for which one can study thoroughly the dispersion curves at the bottom 
of the spectrum, was introduced in uniini. It consists of operators that 
allow for what the authors of usiini call regular factorization: 

(5.4) H = J{^b{DyG{x)b{D)f{x). 

Here b{D) = ^ bjD^^ is a constant coefficients hrst-order linear dif¬ 
ferential operator, the periodic function f{x) is such that f,f~^ G 
and the periodic matrix-function G{x) satishes cqI < G{x) < 
Cil for some 0 < cq < Ci. Then, it was shown in nTi that the bottom 
of the spectrum of if is 0 and it is attained at the quasimomentum 
f = 0, where the band function has a simple nondegenerate minimum. 
As in the previous subsections, this implies 

Theorem 25. (1) The Liouville theorem holds for the operator H 

ig. 

(2) dim Vat (if) = hn,N- 

(3) Every solution u G Yn{H) can be represented as 

u{x) = x^pj{x) 

\3\<N 

with G-periodic functions pj. 
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5.4. Nonselfadjoint second-order operators. Throughout this sub¬ 
section, we again consider the flat case X = M"" only. However, the 
resnlt still holds on general abelian coverings, which does not reqnire 
any signiflcant change in the proof. 

In all the examples that were discnssed so far in this section, we 
considered only selfadjoint operators. There is, however, a class of 
(in general nonselfadjoint) second-order operators in M", which can be 
studied thoroughly and which plays an important role in probability 
theory. 

Consider second-order nniformly elliptic operators on of the fol¬ 
lowing form (EZZl): 

n n 

(5.5) L = — aij{x)didj + bi{x)di + c{x) x G M"" 

i,j=l i=l 

with coefficients that are real and periodic. 

For an operator of this type, an important fnnction A(,^) : M” —^ M 
can be introdnced^®, whose properties were stndied in detail in El El 
[70]. It is defined by the condition that the equation 

Lu = A(^)m 

has a positive Bloch solntion of the form 

(5.6) u^{x) = e^'^p^{x), 

where p^{x) is G-periodic. 

We also need to define the following number: 

(5.7) Ao = maxA(f). 

geR" 

It follows from Eiini that Ao can also be described as follows: 

(5.8) Ao = snp{A G M | 3m > 0 such that {L — X)u = 0 in R""}. 

In the selfadjoint case, Aq is the bottom of the spectrum of the operator 
L. The common name for Aq is the generalized principal eigenvalue of 
the operator L in R"". 

We assemble the information we need abont the fnnction A(^) and 
the generalized principal eigenvalue Aq in the following lemma. The 
reader can find proofs of these statements in EUSHIl and more de¬ 
tailed references in j58j . 

Lemma 26. (1) The value A{^) is uniguely determined for any ^ G 



^®The notation A has been used before in this paper to denote a somewhat 
different -albeit related - object. However, these two notations do not collide in 
this subsection, so the reader should not get confused. 
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(2) The function A{^) is bounded from above, strictly concave, an¬ 
alytic, and has a nonzero gradient at all points except at its 
maximum point. 

(3) Consider the operator 

L(0 = = L{x, D - i^ 

on the torus Then A(,^) is the principal eigenvalue of L{ff) 
with a positive eigenfunction p^. Moreover, A(^) is algebraically 
simple. 

(4) The Hessian o/A(^) is nondegenerate at all points. 

(5) Aq > 0 if and only if the operator L admits a positive (super-) 
solution. This condition is satisfied in particular when c(x) > 0. 

(6) Aq > 0 if and only if the operator L admits a positive solution 
of the form 45.61) . 

(7) Aq = 0 )/ and only if the eguation Lu = admits exactly one 
normalized positive solution in M". 

(8) If c{x) = 0, then Aq = 0 if and only if J h{x)'i({x) da; = 0, 

where i) is the principal eigenfunction of L* on and b{x) = 
{bi{x),... ,bn(x)). In particular, divergence form operators sat¬ 
isfy this condition. 

(9) Let ^ G M”, and assume that u^{x) = e^'^p^{x) and u*_^ are 
positive Bloch solutions of the eguations Lm = 0 and L*u = 0, 
respectively. Denote by the periodic function u^u*_^ . Consider 
the function 

n 

bi{x) = bi{x) -2^aij{x) + [p^{x)]~^ djp^{x)} . 

i=i 

Then Aq = 0 if and only if 

J bi{x)'fi{x) da; = 0 a = 1,..., n. 

Now one can describe Liouville type theorems for an operator of the 
form (ESD assuming that Aq > 0. This assumption implies that the 
operator admits a positive supersolution. 

It is shown in ra that if A(0) > 0, then the Fermi surface can 
touch the real space only at the origin (modulo the reciprocal lattice 
G* = (27rZ)”) and in this case A(0) = 0. 

In fact, we have the following stronger result which extends Lemma 
15 in to nonreal A. In addition, the proof of the statement below 
is more elementary than in that lemma. 
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Lemma 27. Let k = j G C"". If Re \ < A(^), then {k,\) does 
not belong to the Bloch variety of the operator L. Moreover, if 
Re \ = A(^), then {k,X) G B^ if and only if ■j G G*, ^ belongs to the 
zero level set S of A{^), and ImX = 0. 

Proof: Without loss of generality, we may assume that ^ = 0, A(^) = 
0, and LI = 0. Assume that {k, X) G Bl, and let u{x) = e^'^'^p{x) be a 
Floquet solution with a quasimomentum k of the equation 

Lu = Xu in M”, 

where Re A < 0. Take the complex conjugate 

Lu = Xu in M”, 

and compute 

n n 

L(|m|^) = uLu + uLu — 2 aijUxfixj = 2i?e A|n|^ — 2 aijUxfixj- 

i,j=l i,j=l 

Notice that for each G C", we have 

n n n 

aijCiCj = OijReCiReCj + aijImC,iImC,j > 0. 

i,j=l i,j=l ij=l 

Therefore, 

(5.9) L{\u\‘^)<0. 

Thus, Inp = |p(a;)p is a periodic nonnegative subsolution of the equa¬ 
tion Lu = 0 in R”'. By the strong maximum principle \u{x)\^ = 
constant. In particular, L{\u\‘^) = 0. Since we have equality in (15.9j) . it 
follows that i?e A = 0 and 

n n n 

^ ^ ^ijUx^Uxj ^ ^ aijRe Ux^Re Uxj T ^ ^ aijlmux^lrnuxj 0. 

i,j=l i,j=l i,j=l 

It follows that u = constant and therefore, is a periodic function. 
Consequently, 7 G G*. Since LI = 0, it follows that ImX = 0. □ 

Lemma 1771 and Theorem ITHI imolv that if A(0) > 0, then the Liouville 
theorem holds vacuously, and the equation Ln = 0 does not admit any 
nontrivial polynomially growing solution. 

Suppose now that A(0) = 0 (i.e., 0 G Fl,r)- Then by Lemma 1771 
Fl,r = {0}. Moreover, Lemma 1^ implies that if Aq > 0, then the point 
A: = 0 is a noncritical point of the dispersion relation, and if Aq = 0 
then /c = 0 is a nondegenerate extremum. In terms of Theorem we 
have in the first case /q = 1, while in the second Iq = 2. Moreover, in 
both cases r = 1. Theorems ITHl and ITHl imolv now that the Liouville 
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theorem holds, and every solution u G Yn{L) is representable in the 
form (ca). The dimension of the space Yi^{L) is equal to hn,N in the 
case when Aq = 0, and to qn-i,N when Aq > 0. 

We summarize these results in the following statement: 


Theorem 28. Let L be a periodic operator of the form such that 
Aq > 0. Then 

(1) The Liouville theorem holds vacuously if A{0) > 0, i.e., the 
equation Lu = 0 does not admit any nontrivial polynomially 
growing solution. 

(2) //A(0) = 0 and Aq > 0, then the Liouville theorem holds for L, 

dimVAr(L) = qn-i,N = 

and every solution u G Yp^{L) of the equation Lu = 0 can he 
represented as 

u{x) = x^pj{x) 

\'j\<N 


( n + N — 1 '\ 

V ^ ) 


with G-periodic functions pj. 

(3) //A(0) = 0 and Aq = 0, then the Liouville theorem holds for L, 


dim VAr(L) = hn^jsf = 


n + N 
N 


n + N — 2 \ 
N-2 ) ’ 


and every solution u G Yn{L) of the equation Lu = 0 can he 
represented as 


u 


(x) = ^ x^pj{x) 
\j\<N 


with G-periodic functions pj. 


It is interesting to notice that for selfadjoint operators, one never 
faces the case / = 1 at the spectral edges, as it does happen above 
when A(0) = 0 and Aq > 0. 


6. Liouville theorems for overdetermined elliptic 

PERIODIC SYSTEMS AND FOR HOLOMORPHIC FUNCTIONS 

We will show here how the techniques and results of the preceding 
sections can be applied to some overdetermined elliptic systems, in par¬ 
ticular to the Cauchy-Riemann operators on abelian covers of compact 
complex manifolds. The main construction stays the same, so we will 
be brief in its description. 
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We consider as before an abelian covering X i—> M of a compact 
Riemannian manifold. Let also Lj (j = 0,1,...) be finite dimensional 


smooth vector bundles on M equipped with a Hermitian metric and 
Cj be the sheaves of their smooth sections. Suppose that we have an 
elliptic complex of differential operators 


Pi 


Co ^ ClC 


P2 


We can lift this complex to an elliptic complex on X for which we will 
use the same notations: 


( 6 . 1 ) 


Co^Ci 


Pi 


C, 


P2 


Notice that the lifted bundles have a natural G-action with respect to 
which the operators are periodic. 

We will be interested in the spaces VAr(P) of the global solutions 
u{x) of the equation Pu = 0 on X that have polynomial growth in the 
sense that for any compact X G X 

(6.2) \\u\\L'2{gK,Lo) ^ + \g\)^ WgeG 

for some C and N. As before, without loss of generality we can assume 
that G = Z”. The Floquet-Gelfand transform reduces the elliptic 
complex (EH) and its dual to the direct integrals with respect to the 
quasimomentum k of elliptic complexes on M 

(6.3) 0 - ..., 


and 

(6.4) 


0 ^ r- 


where are appropriately dehned analytic Banach vector bundles 

over C” (a more detailed setup of the relevant spaces can be found 
for instance in [H].). As previously, solutions of polynomial growth of 
Pu = 0 after applying the Floquet-Gelfand transform become func¬ 
tionals on G°°(T”,T^) orthogonal to the range of the operator 


(6.6) 


C“(T", J^^) '"'S^ C“(T", J^"). 


Notice that due to the ellipticity of the complex and the compactness of 
the base, it is Fredholm as a complex of operators between the appro¬ 
priate Sobolev spaces (e.g., im Vol. Ill, Sections 19.4 and 19.5], jTHl 
Section 3.2., in particular Theorem 13], |841 Section IV.5]). Hence, the 
operators P*{k) on M belong to the set of the right semi-Fredholm 
operators between the corresponding Sobolev spaces of sections. This 
means that they have closed ranges of hnite codimension, while having 
inhnite dimensional kernels. It is also clear that as in the elliptic case 


LIOUVILLE THEOREMS ON ABELIAN COVERINGS 


39 


before, P{k) depends analytically (in fact, polynomially) on k. The 
Fermi surface Fp of P is introduced as the set of all values of k for 
which the equation P{k)u = 0 has a nonzero solution^^. As before, it 
coincides with the set where the dual operator P*{k) has a nonzero cok¬ 
ernel. This enables us to carry over the considerations of the previous 
section to prove the following result: 

Theorem 29. The following conditions are equivalent: 

(1) The (real) Fermi surface Fp^r is finite. 

(2) For some N >0 the dimension dN of the space VAr(P) of solu¬ 
tions of polynomial growth of order N is finite (i.e., the Liouville 
theorem of order N holds). 

(3) The dimension of Yp[{P) is finite for any N > 0 (i.e., the 
Liouville theorem holds) and < d^qn^N- 

If one of the above conditions holds, then any solution of polynomial 
growth is a linear combination of Floquet solutions. 

Proof: Since any point k G Fp,R provides a bounded Bloch solution 
and those for different values of k (modulo the dual lattice) are linearly 
independent, the implications (3) ^ (2) ^ (1) are obvious. Let us now 
assume (1) and prove (3). According to Lemma El all functionals of 
order N orthogonal to the range of (p31) can be expressed as linear com¬ 
binations over k G Pp^R of the form (14.411 . The same Lemma together 
with Lemma ITUl and the hniteness of Pp^R show the hnite dimensional¬ 
ity of this space and the Floquet representation for all solutions of this 
type. □ 

We apply the above theorem to the following interesting case. 

Theorem 30. Let X ^ M be an abelian cover of a compact (complex) 
analytic manifold M. Let X be equipped with the lift of an arbitrary 
Riemannian metric from M and p be the corresponding distance. For 
N = 0,1,..., let Ais[{X) be the space of all analytic functions on X 
growing not faster than C(1 -|- p{x))^. Then: 

(1) The real Fermi surface of the d-operator defined on functions 
on X contains only the origin: P^r = {0}. 

(2) Aq{X) consists of constants. 

(3) For any A^ > 0, dimAAr(X) < q^^pi < oo. All elements of 
the space Ajv(W) are holomorphic Floquet-Bloch functions with 
quasimomentum fc = 0. 


^^It can be shown similarly to the elliptic case that the Fermi surface, as in the 
elliptic case, is an analytic subset of C”. 
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(4) There exists a finite family of holomorphic functions {fj} on X 
such that all elements of An{X) are polynomials of {fj}. 

Proof: Let k G be in Fg^. This means that there exists a nonzero 
analytic function u{x) of Bloch type, i.e., automorphic with a unitary 
character 7 ^. of the group G. The absolute value of u{x) is then G- 
periodic and thus can be pushed down to M. Due to the compactness 
of M, it must attain its maximum. Hence, the function u{x) attains 
somewhere on X the maximum of its absolute value. The standard 
maximum principle now implies that u = constant. This means in 
particular that k = 0. Hence, (1) holds. 

Representation (D of Lemma fT?)l implies now that Aq{X) consists 
of constants. Indeed, this space is comprised of analytic functions on X 
automorphic with respect to unitary characters of G. As we have just 
shown, every such automorphic function is automatically constant. 

Let us extend the d operator on X to the elliptic Dolbeault complex 
(e.g., IHH Chapter IV, Examples 2.6 and 5.5]): 

( 6 . 6 ) 

where A^’’^ is the sheaf of germs of smooth (p, q')-forms on X. The 
conditions of Theorem 1221 are now satished, which proves the statement 
(3). 

Finally, the implication (3)^(4) was established in ^2]- D 

7. Liouville theorems on combinatorial and quantum 

GRAPHS 

Liouville type theorems have also been studied intensively on dis¬ 
crete objects (graphs, discrete groups) (e.g., jlSlISHI)- Our technique is 
insensitive to the local structure of the base of the covering, and hence 
it can also be applied to abelian co-compact coverings of graphs. 

Let T be a countable noncompact graph with a free co-compact ac¬ 
tion of an abelian group G of a free rank n. Denote by T the hnite 
graph T/G. Consider any periodic hnite order difference operator P 
on T. One can define all the basic notions: solutions of polynomial 
growth, Floquet solutions, quasimomenta, Fermi variety, etc., exactly 
the way they were defined before. The same procedure as before leads 
to the following result that we state without a proof: 

Theorem 31. (1) If the Liouville theorem of order N > 0 for the 

equation Pu = 0 on T holds true, then it holds to any order. 

(2) In order for the Liouville theorem to hold, it is necessary and 
sufficient that the real Fermi surface (i.e., the set of unitary 
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characters x for which the equation Pu = 0 has a nonzero x- 
automorphic solution) is finite. 

(3) If the Liouville theorem holds true, then: 

(a) Each solution u G Vn^P) can be represented as a finite 
sum of Floquet solutions of order up to N. 

(b) If dN = dimY n{,P) , then for all N eN, we have 


dN < < oo , 


where qn,N is the dimension of the space of all polynomials 
of degree at most N in n variables. 

(c) Assume that for each real quasimomentum k G the 
conditions of Theorem^^ are satisfied. Then for each N > 
0 the dimension dpf of the space Yn{P) is equal to 




n + N — lo{k) 
N-lo{k) 


n + N 
N 


Here r^ and Ifik) are respectively the multiplicity of the 
zero eigenvalue and the order of the first nonzero Taylor 
term of the dispersion relation at the point k G -Fp^R. The 
terms in this sum are the dimensions of the spaces of 
harmonic polynomials, and polynomially growing solutions 
can be described in terms of these polynomials analogously 
to Theorem m 


The same can also be done for the so called quantum graphs, where 
a differential (rather than difference) eqnation is considered on a graph 
that is treated as a one-dimensional singnlar variety rather than a 
pnrely combinatorial object (see the details in in2]-insi)- A qnantnm 
graph r is a graph with two additional strnctnres. First of all, the graph 
mnst be metric, i.e., each edge e mnst be snpplied with a (hnite in onr 
case) length /e > 0 and correspondingly with an ‘arc length’ coordinate. 
This allows one to do differentiations and integrations and to dehne 
differential operators on T. The simplest ones are —dfix/dx^ + V{x), 
where differentiation is done with respect to the edge coordinates, and 
V{x) is a snfficiently nice (e.g., measnrable and bonnded) potential. In 
order to dehne snch an operator, one needs to impose some bonndary 
conditions at vertices. All snch ‘nice’ conditions have been described 
(e.g., 1211 EHES] )• The simplest example is the so called Nenmann (or 
Kirchhoff) condition that reqnires that the fnnctions are continnous at 
the vertices and that the snm of ontgoing derivatives at each vertex 
is eqnal to zero. A quantum graph is a metric graph eqnipped with a 
selfadjoint operator P of the described kind. 
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Let now L be a noncompact quantum graph with a free co-compact 
isometric action of an abelian group G of a free rank n, and let L be 
the compact quantum graph T/G. We assume that the Hamiltonian 
P of the quantum graph is G-periodic. The basic notions (solutions 
of polynomial growth, Floquet solutions, quasimomenta, Fermi variety, 
etc.) can again be defined for this situation. In particular, we consider 
the operators P{k), and assume that for each k E the spectrum of 
P{k) is discrete. 

Theorem 32. The statements of Theorem\^ hold for the quantum 
graph F under the conditions just described. 

8 . Further remarks 

(1) In most parts of the paper we assumed that all the coefficients 
of the operators P and P* are G“-smooth. This assumption 
was made for simplicity only, and in fact one does not need 
such a strong restriction (e.g., see the discussion in m Remark 
6.1] and Section 3.4.D]). For example, for the results of 
Theorem|2Hlto hold, it is sufficient (but not truly necessary) that 
the coefficients of L and L* are Holder continuous. It is clear 
that the conditions on the coefficients could be significantly 
relaxed if the operators were considered in the weak sense, or 
by means of their quadratic forms. This does not change the 
general techniques of the proofs, which rely on the analytic 
dependence with respect to the quasimomenta, rather than on 
the particular regularity properties of the coefficients. One only 
needs to guarantee the compactness of the resolvents of the 
operators P{x)- More precisely, we need that both the operator 
Pm and its dual Pfj dehne Fredholm mappings between the 
Sobolev space H^{M) and L 2 {M), and this condition can be 
weakened further. 

(2) A stronger statement than Theorem El holds: the existence of 
a subexponentially growing solution also implies the existence 
of a Bloch one. This is an analog of the Shnol’s theorem (see 
IZ3ISZIIZ3). It is provided in Theorem 4.3.1 of ISH for the case 
of periodic equations in M”, however carrying it over to the 
case of more general abelian coverings does not present any dif¬ 
ficulty. An analogous theorem holds for periodic combinatorial 
and quantum graphs |HH] . 

(3) The interesting feature of the main results of the present paper 
is that they relate the Liouville property to the local behavior of 
the dispersion relations near the ‘edges’ of the spectrum of the 
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operator. This threshold behavior is responsible for many phe¬ 
nomena, in particular homogenization jsiiiniiisiiinii23i2niii2i, 
structure of the impurity spectrum arising when a periodic 
medium is locally perturbed 11211111, Anderson localization, and 
others. There are less explored issues of this kind, for instance 
the behavior of the Green’s function |E1 Ea EDI EH, integral 
representations of solutions of different classes, and finally Li- 
ouville theorems. The last two were also looked upon in our 
previous paper j2H|- The results of the papers quoted above, 
as well as of this paper, deduce properties of solutions from 
an assumed spectral edge behavior of the dispersion relation. 
However, there are only extremely few results concerning such 
precise spectral behavior for specific (or even generic) opera¬ 
tors. The notable exception is the bottom of the spectrum, 
where much is known (e.g., EHIZHIIHI and recent advances in 
Ha im EO])- An initial study of the internal edges was con¬ 
ducted in the interesting paper gni, however the knowledge in 
this case remains far from being satisfactory. It is also appro¬ 
priate to mention that the study of analytic properties of the 
Fermi and Bloch varieties, even for ‘generic’ operators, happens 
to be a very hard problem (e.g., |2aEnillZl), and there are sev¬ 
eral unproven conjectures about the generic behavior of these 
varieties (see for example, 01711). 

(4) As it was mentioned above, the behavior at the bottom of the 
spectrum is responsible for homogenization. Our results on Li- 
ouville theorems indicate that in some sense some kind of “ho¬ 
mogenization” exists also at the interior edges of the spectrum. 
A recent progress has been made towards precise understanding 
of the meaning of this ‘homogenization’ in m for ODEs and 
[TT)] for PDEs. 

(5) An unexpected and counter-intuitive result of our study of Li- 
ouville theorems is that they depend on the order Zq of the 
leading term Xi^ of the dispersion relation only, rather than on 
its full structure. For instance, if Iq = 2, the dimensions are the 
same as for the Laplace operator, even if A 2 , the Hessian of the 
dispersion relation, is degenerate (e.g., the case of the operator 
Mt with t = 2a/ 3 in Sectional). Normally in homogenization 
theory, the whole term is important. 

(6) Although our results are complete for the case of non-multiple 
spectral edge (r = 1 in terms of Theorem fTT^ and II bj) . in the 
multiplicity case when r > 1 we require that the determinant of 
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the leading term A/g of the Taylor expansion is not identically 
equal to zero. Simple examples show that without this non¬ 
degeneracy condition, the proof of M.fij) (and probably the result 
as well) does not work anymore. Indeed, suppose that we are 
dealing with a single point k of the Fermi surface and that the 
matrix X{k) is diagonal, with each diagonal entry Xjj{k) having 
its own leading term Xi {k). In this case, it is not hard to check 
that the summand in (HU) needs to be replaced by 


( 8 . 1 ) 


^ / n-\- N \ f n ^ N — Ij \ 

vK ^ ; V N-i, )_ 


in order for the result to hold. One can also see that in the 
non-selfadjoint case the non-degeneracy condition in its present 
form essentially forces the algebraic and geometric multiplici¬ 
ties of the zero eigenvalue to coincide. These dehciencies can be 
mollihed by allowing more general ways of choosing the lead¬ 
ing term A/^, e.g., in the ‘Douglis-Nirenberg sense’ analogous 
to the corresponding dehnition of ellipticity j2H|- However, at 
the present moment it is not clear to the authors how to avoid 
having any non-degeneracy conditions in the case of a multi¬ 
ple eigenvalue. Fortunately, in most known to us cases that 
are interesting for applications, one expects multiplicity to be 
absent. 

(7) The result of Theorem 1201 for the case of Kahler manifold can be 
extracted from the elliptic case, by switching the consideration 
to harmonic functions. Such result was obtained in ^O]- This, 
however, does not work in general and the technique presented 
in this text is needed. 

(8) It is interesting to mention that the dimension (and in fact the 
geometry) of the underlying manifold X does not explicitly en¬ 
ter in our main results. It dehnitely influences the geometry of 
the dispersion curves and hence the Liouville theorems. How¬ 
ever, on the surface it looks like the group G = Z"' matters more 
than the manifold. In particular, one can easily cook up a two 
dimensional covering X —> M with an arbitrarily large n (take 
for example the standard 2-dimensional jungle gym JG^ in M”, 
see I7SI). Then in the results that concern the dimensions of 
the spaces of polynomially growing solutions, one sees mostly 
the influence of n. In particular, one can get large dimensions 
of these spaces in such a manner. The reason is that, as it has 




LIOUVILLE THEOREMS ON ABELIAN COVERINGS 


45 


been mentioned before, onr Liouville theorems are of ‘homoge¬ 
nized’ nature, and being seen from afar, the covering manifold 
X looks like [SllEH]- 
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